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摘 要 : 为 了 改善 瑞 利 衰落 信道 单 径 自 相关 统计 特性 ， 提 出 了 一 种 改进 的 多 径 非 相关 无 线 信 道 仿真 模型 。 该 模型 首次 
在 荣 特 卡 罗 法 (MCM) 的 到 达 角 (AOA) 中 引入 多 径 因素 k, 改进 了 模型 AOA 的 取 值 方法 。 通 过 对 信道 模型 参数 AOA、 
k、 低 频 振荡 器 个 数 N 以 及 均 方 误差 (MSE) 的 仿真 分 析 ， 结 果 表 明 ， 该 模型 在 相关 统计 特性 方面 与 参考 模型 相 比 具 
有 极 好 的 吻合 性 ; 其 次 ， 与 MCM 相 比 ， 在 实现 开销 相同 的 情况 下 ， 改 进 的 新 模型 在 最 大 允许 时 间 间 隔 内 ， 衰 落 信道 
单 径 自 相 关子 数 的 MSE 约 为 MCM 的 四 分 之 一 ， 精 度 有 了 明显 的 提高 ， 能 更 加 准确 地 描述 无 线 信 道 。 新 模型 为 无 线 信 
道 仿 真 器 的 设计 、 频 率 选 信道 的 建 模 等 提供 了 一 种 新 的 方法 。 
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wireless channel model 
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Abstract: In order to improve the statistical properties of multipath fading and more accurate description of wireless channel. 
This paper proposed an improved simulation model based on Monte Carlo method (MCM) for uncorrelated multiple Rayleigh 


fading channels. Based on the angle of arrival (AOA) , the multipath element (k) is introduced into the new simulation model, 


加 入 which improved the method of evaluating AOA. Through the simulation analysis of the channel model parameters AOA, k, the 


number of low frequency oscillator and the mean square error (MSE) , the results show that a pretty good agreement statistical 


property of the new model and those of the underlying reference model has been observed. Compared with MCM, the new 
model has a significant improvement in the autocorrelation properties in the same conditions and the MSE is about a quarter of 
the original. The new method is useful for the design of wireless channel simulator, frequency-selective channel. 
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复 高 斯 随机 过 程 ， 而 基于 SOS 的 信道 模拟 方法 能 够 模拟 瑞 利 、 
莱 斯 等 多 种 无 线 衰落 信道 , 可 以 根据 需要 产生 特定 长 度 的 波形 ， 
随 着 多 输入 多 输出 (multiple input multiple outpub MIMO) 、 且 具 有 易于 实现 和 精确 有 效 的 优点 ， 在 实际 仿真 模拟 中 应 用 更 
分 集 技术 、 宽 带 无 线 网 络 通信 系统 的 发 展 与 应 用 ， 无 线 信道 模 。 ”为 广泛 。 基 于 SOS，Jakes 在 Clarke 提出 经 典 信道 参考 模型 后 ， 
拟 在 现代 无 线 通 信 设 备 研发 与 测试 的 过 程 中 所 占 的 比重 越 来 越 ”提出 的 简化 衰落 模型 在 信道 平稳 方面 以 及 统计 独立 性 方面 存在 
大 。 其 中 ， 模 拟 数 字 通 信 的 多 径 衰 落 信道 时 ， 由 于 真实 的 无 线 。 着 不 足 之 处 ， 产 生 的 信号 也 不 是 各 态 历 经 的 ， 各 径 信 道 之 间 相 


通信 信道 复杂 多 变 ， 能 否 准 确 地 描述 无 线 信 道 直 接 决 定 着 无 线 。 关 性 强 。 

通信 质量 的 好 坏 ， 因 此 ， 生 成 多 个 单 径 自 相关 统计 性 能 好 且 各 为 了 提高 单 径 信道 的 自 相 关 性 ， 减 小 信道 与 信道 之 间 的 相 

径 不 相关 的 衰落 信道 对 无 线 通信 的 性 能 进行 研究 依然 是 无 线 通关 性 , 于 是 有 了 文献 [5~12] 改 进 的 模型 , 但 是 各 模型 仍 有 很 多 不 

信 信 道 建 模 的 重点 之 一 11。 足 之 处 : 其 中 ，Pop 等 人 在 文献 [5] 中 引入 随机 平均 分 布 相位 ， 
对 于 常见 的 瑞 利 衰落 信道 建 模 ， 分 为 滤波 法 和 正弦 波 羡 加 ”解决 了 信道 非 平稳 性 问题 ， 但 仍 不 能 生成 多 径 非 相关 信道 


法 (sum-of-sinusoids, SOS) 。 滤波 法 运算 量 大 且 需 要 产生 独立 Zheng 等 人 文献 [6,7] 中 提出 的 不 确定 性 模型 ,大 大 改善 了 瑞 利 误 
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落 信 道 正 交 分 量 间 的 互相 关 性 ， 并 且 Zheng 模型 
计 特 性 ,适用 于 多 径 瑞 利 衰落 信道 ,但 是 Zheng 模型 的 到 达 角 只 
考虑 了 同 径 信 道 的 随机 分 布 特性 ， 并 没有 包含 多 径 因 素 ; Li 等 
人 在 文献 [8] 中 合理 地 设计 多 普 勒 频 移 ,实现 了 多 径 独 立 瑞 利 衰 
落 信道 的 仿真 ， 但 到 达 角 不 是 随机 分 布 ， 在 文献 [9] 中 MCM 模 
型 虽然 到 达 角 是 随机 分 布 ， 但 是 相 比 之 下 ， 人 
能 差 。 而 文献 [10~12] 改 进 的 模型 都 是 在 等 距 法 (MED )、 
只 法 (MEA ) 或 者 是 精确 多 普 勒 扩展 法 (MEDS ) 基础 上 引入 
多 径 因 素 提出 的 ， 且 模型 的 自 相 关 ， 互 相关 等 相关 统计 特性 与 
参考 模型 仍然 存在 着 很 大 的 差别 。 
综 上 所 述 ， 针 对 多 径 非 相关 无 线 信道 的 特点 ， 分 析 了 文献 


有 很 好 的 统 
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式 (6) 中 Cix(t)，Czx(t) 分 别 表 示 第 大 径 的 同 向 分 量 和 正 
交 分 量 , k=1,2,…,M ，M 为 不 相关 波形 的 个 数 ; 式 (7) 中 ce、 
finx 和 i 分别 表示 第 径 各 支 路 的 多 普 勒 系数 、 多 普 勒 频 
率 和 随机 相位 。c,， N =N ,n=1,2,.…,N,(i=1,2), 
Ji= 太 sin ，0sV(02rl ， 大 表示 第 大 径 。 
传播 理论 可 知 ， 接 收 信号 是 由 所 有 散射 波 的 矢量 县 加 而 
ee 影响 信道 模型 的 统计 特性 。 
文献 [13] 根 据 cs、 en 
型 分 为 四 大 类 ， 和 (cn 是 确定 变量 ，fi, 和 ,是 
随机 变量 ) 是 高 效 的 瑞 利 信道 仿真 模型 。 类 型 II 中 有 两 种 典型 
的 仿真 模型 : MCM 模型 和 Zheng 模型 。MCM 模型 的 统计 特性 


=\J2ci/N ， 


中 己 有 的 衰落 信道 模型 与 各 个 模型 不 足 之 处 的 基础 之 上 ， 本 文 
首次 尝试 在 传统 MCM 模型 的 AOA 中 引入 多 径 因 素 ， 提 出 了 
一 种 新 的 AOA 模型 ,改进 了 MCM 的 实现 方法 。 本 文 提 出 的 改 
进 方法 在 各 径 的 自 相 关 特 性 上 得 到 了 明显 提高 ,而 且 多 径 信道 
之 间 不 相关 ， 能 更 加 准确 地 描述 无 线 衰落 信道 。 


1 。 瑞 利 衰落 信道 参考 模型 


基于 SOS 的 瑞 利 衰落 信道 仿真 模型 ， 
弦 波 短 加 实现 的 “ 


是 通过 有 限 数 
模型 如 式 (1) 和 (2) 所 示 。 


的 正 


cO=caD+jcz) (1) 
ci(t) = yes cos(2nf. ,1 + 0.,) O) 


式 (D) 中 ct(D) ， ca(D 分 别 为 同 向 分 量 和 正 交 分 量 ， 式 (2) 
中 i=1,2 ，N 为 模拟 瑞 利 信道 需要 的 低频 振荡 器 个 数 , 参数 fi,、 
0., 分 别 为 各 支 路 离散 多 普 勒 频率 和 多 普 勒 相位 ， 参 数 c., 为 信 
道 多 普 勒 系数 ， 如 式 (3) 和 式 (4) 所 示 。 
ca = ci/ G) 
fin = fa sin(a;,,) (4) 
式 (4) 中 为 最 大 多 普 勒 频 移 ， 4 为 各 支 路 信号 的 AOA。 
该 模型 对 应 的 二 阶 统计 参量 表示 为 
Ruan(r)=Reo()=aJo(wz) 
Rae,(T)= Roalt)=0 (5) 
Ree(T) =20°, (wt) 
式 (5) 中 由 表示 零 阶 第 一 类 Bessel 函数 ，o? 为 瑞 利 衰落 


的 方差 ，R6。 (7) 表示 单 径 信道 的 同 向 分 量 自 相 关 函 数 ， Re (7) 
表示 单 径 信道 的 正 交 分 量 互相 关 函 数 ; Re(z) 表示 单 径 信道 的 
自 相 关 函 数 。 

2 ”改进 的 MCM 模型 

2.1 改进 MCM 模型 的 AOA 表示 

根据 瑞 利 衰落 信道 的 参考 模型 表示 方法 ， 在 AOA 中 引入 


多 径 因 素 后 , 非 相关 多 径 瑞 利 衰落 信道 
Cue(D=Cu(D+jCze(D (6) 


N, 
Cix(D) = Dc Cos(2f st + Oi) (7) 
n=] 


模型 表示 为 式 (6) 和 (7) 。 


与 参考 模型 的 统计 特性 偏差 较 大 ; Zheng 模型 相对 而 言 统 计 特 
性 与 参考 模型 更 接近 ， 但 不 具备 各 态 历经 性 。 且 各 模型 都 不 有 具 
有 多 径 因素 ， 具 体 AOA 如 式 (8) 和 (9) 所 示 。 


MCM 模型 : an = 也 的 (8) 
2xn—x+g 
Zh 弄 J. a 9 
eng 模型 RE 0) 


式 (8) (9) 中 n=1,2,…， 
peU(-A,A]。 
式 (9) 可 知 ， 乡 取 不 同 值 时 ，Zheng 模型 的 所 有 AOA 取 


值 在 rr-DAmgA 


Ni=L2) ; wu,eU(0,]] 


随机 分 布 , 当 n=N 时 , AOA 有 最 大 值 


zj2 ,如 图 (1) 所 示 。x 轴 取 值 [1,N] , 且 为 整数 ,了 轴 取 值 (0,2] 。 


图 1 AOA 随机 取 值 区 域 


分 析 比 较 上 述 式 (8) 与 (9) 的 AOA 取 值 范围 与 特点 ， 在 
MCM 模型 的 AOA 中 引入 多 径 因 素 重新 对 AOA 取 值 ， 使 
之 : 


内 随机 分 


T 
tr 


3 无论 k 取 何 值 ，AOA 在 70" ,m7w4wl 范 上 


布 ， 进 而 改善 MCM 信道 模型 的 统计 特性 ; 
b) 当 n=1 时 , AOA 大 于 0, n=NN 时 , AOA 小 于 等 于 7/2; 
co) 改变 k 的 值 可 以 快速 生成 所 需 的 不 相关 瑞 利 波形 。 
新 模型 的 AOA 如 式 (10) 所 示 。 可 以 看 出 各 模型 的 根本 
别 在 于 AOA 的 取 值 方法 不 同 。 
新 模型 : 


[xl 
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工 ， Ax(n—1) 
2 ”2(V-D 


(10) 


-|o- M0 


5 
式 (10) 中 wsv[0m2) i=1,2; wu, EeU(0,l]; k=1,2,.…,M ; 

M 为 生成 的 不 相关 波形 的 个 数 。 

2.2 改进 的 MCM 模型 分 析 
根据 参考 模型 ， 该 模型 的 第 

相关 函数 ACF 如 式 (11) 所 示 " 


k 径 瑞 利 衰落 波形 Cix(D 的 自 


lim im 元 六 Ci (DCi (t+ Td 


= E[Cix(1)Cix(t +7)] 
N 
im Dcos(2nf sat 二 Ok ): 
en 


N 
cin 2 cos (Qnfiss (t+ 7) + Oa)] 
n=] 


Doosan) (11) 
式 (11) 中 i=1,2,k=1,2,…,M ;7 为 延迟 时 间 ; M 为 生成 
的 不 相关 波形 的 个 数 。 
同 理 , 瑞 利 衰落 波形 Cix(7) 第 k 径 的 互相 关 函 数 CCF 如 式 
(12) 所 示 
reic, = 
lim 地 | Cix(f)Cju(t +T)dt 
= E[Cix(f)Cju(t +7)] 
=0 (eh ps) (12) 
即 当 iz# jG j=1,2) 和 klk,l=1,2." M), Cix(t), 


普 勒 频率 无 重 值 ，Cix (71) ， 


Cz.1(?) 的 离散 多 CD 的 互相 关 函 


数 CCF 值 为 零 。 
为 了 满足 式 (12) ， 模 型 各 径 信 道 间 的 CCF 期 望 值 为 零 ， 
AOA 应 满足 
cos(a, 1) +sin(a,,) 
os *) cos(ans) [kh (13) 


sh 


式 (13) 中 k,m=1,2,.…,M ， M 为 不 相关 波形 的 个 数 ， 
n,1=1,2,…,N ，N 表示 需要 的 低频 振荡 器 个 数 。 
改进 模型 的 AOA 满足 上 面 的 约束 条 件 ， 将 上 式 (10) 代入 


到 式 (11) 中 ， 改 进 模型 的 自 相关 函数 ACF 见 式 (14) 。 
工 和 (1 一 ]) 
Feos 27nf, sin < 2 (14) 
-|e- 10 -+ 3h, 

3 “仿真 分 析 
3.1 NN 值 影响 仿真 分 析 

以 改进 的 MCM 模型 、MCM 模型 与 参考 模型 的 单 径 信道 
同 向 分 量 自 相关 函数 ACF 的 MSE 作为 判断 标准 , 分 析 N 的 不 


可 


取 值 对 模型 精度 的 影响 ， 即 
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(15) 


由 2 ] 
MSE=E IE 一 Re 


在 入 取 不 同 值 时 ， 改 进 MCM 与 MCM 的 i MSE 值 如 图 
2 所 示 。 从 图 中 可 以 很 直观 地 看 出 ， 与 MCM 相 比 ， 无 论 入 取 
何 值 ( 如 N=11 时 ，MCM 的 值 为 0.045 39， 而 改进 MCM 的 
值 只 有 0.01031)， 改 进 MCM 得 到 的 结果 都 明显 要 好 ， 而 且 新 


模型 相 较 于 原先 的 模型 ， 相 关 性 质 收敛 快 。 再 一 次 证 明了 新 模 
型 能 更 加 准确 地 描述 无 线 衰 落 信道 
0.1F 
0.08 - 
Al ee 
器 006- 0 
a 
0.04 - 
0.02 | 
0 . 
0 +10 


图 2 不 同 NM 取 值 的 均 方 误差 


对 于 不 同 的 N 值 ， 改 进 MCM 的 单 径 
值 如 表 1 所 示 。 


相关 函数 的 MSE 


表 1 不 同 入 值 的 MSE 
N 3 9 10 11 12 


MSE 0.02436 0.02255 0.01526 0.01031 0.0090 


随 着 N 值 的 增 大 ， 自 相关 函数 的 MSE 逐渐 减 小 ， 综 合 
虑 精度 、N 值 曲线 的 斜率 变化 和 实现 成 本 ， 本 文选 取 N =11 作 
为 仿真 条 件 。 
3.2 AOA 值 仿 真 分 析 
式 (13) 知 ， 当 AOA 值 wesv[or2) 时， 满足 


Ort a A/2 
Qk 天 dm 1 关 n 或 km (16) 
式 (16) 中 Kk,m=1,2,…,M ， 171 =12,……, N 


当 N=11，M =2 时 ,根据 式 (10) 随机 生成 两 径 信 道 的 
AOA, 分 布 及 正 余弦 函数 值 如 图 3 所 示 。 图 中 横 坐 标 表示 AOA 
值 wx 的 相位 取 值 。 可 以 看 出 各 组 值 均 不 相等 旦 满足 式 (16)。 


1 t 
各 
0.8 
0.6 记 
起 晤 | we # Sinal 
Sin a, 
cos al 
0.2 
cos ao 
0 
0 0.1 0.2 0.3 0.5 
相位 /区 


图 3 两 组 AOA 取 值 及 正 余弦 值 
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3.3” 自 相关 性 能 仿真 分 析 改进 模型 每 一 径 的 ACF 的 精度 都 有 着 巨大 的 提升 .图 6 中 


为 了 直观 上 了 解 改 进 模型 的 性 能 ， 且 有 可 比 性 ， 选 取 x 轴 表 示 k=1,2,…,M 的 各 径 ，y 轴 采 用 fr 作为 标准 量度 ，z 
N=11, o =1, f=100Hz, M=5。. 取 MM 中 的 第 1 径 与 MCM  ， 轴 表 示 各 径 自 相关 ACEF 的 MSE。 
和 参考 模型 进行 仿真 对 比 ， 可 以 发 现 改 进 的 新 模型 的 自 相关 函 表 2 列举 了 不 同上 值 对 应 的 改进 MCM 的 MSE 值 和 MCM 
数 ACF 与 参考 模型 的 ACF 具有 极 好 的 吻合 性 ， 如 图 4 所 示 ， 的 MSE 值 。 分析 表 2 中 的 数据 可 以 看 出 ， 改 进 模型 单 径 信道 
图 中 以 fr 作为 模 坐 标的 标准 量度 ，7 的 最 大 延 时 时 间 取 ”的 每 一 径 的 自 相 关 MSE 的 精度 均 有 极 大 的 提 到 , 约 为 MCM 的 


N/(2f,) 。 四 分 之 一 。 
表 2 不 同 大 值 的 MSE 
| \ erence 径 1 2 3 4 5 
| MCM 
' i 改进 MCM MCM 0.04539 0.04682 0.04835 0.04787 0.04599 
加 改进 MCM 0.01031 0.01062 0.01192 0.01224 0.01213 
| A 
, HH \ A 内 A /~ a Ne 
中 3.5 非 相关 信道 时 域 波形 
1 i 9 [4 4 j | es ee 
oe VY YY 仿真 给 出 了 M =5 实验 中 信道 的 第 1 条 、 第 3 条 和 第 5 条 
时 衰落 信道 的 时 域 包 络 波形 ,如 图 7 所 示 。 从 直观 上 就 可 以 看 出 ， 
i 2 该 三 条 信道 的 包 络 不 具有 相关 性 ， 且 有 较 好 的 非 相 关 瑞 利 衰落 
Ry 信道 的 特点 。 
图 4 各 模型 自 相 关 函 数 
改进 模型 与 MCM 和 参考 模型 的 MSE 如 图 5 所 示 。 从 图 SF 
中 可 以 看 出 ， 在 上 面 的 实现 条 件 下 ， 在 最 大 允许 时 间 间 隔 内 改 of A YN A AM 1/ A 
进 模型 的 误差 仅 为 MCM 的 四 分 之 一 , 性 能 得 到 了 极 大 的 提升 。 Wo A A A 
mo 3h! iN | | | 
0.06 专 妈 -10 | 
| / k=5 
0.04 vf VV 
| 200 0.02 0.04 0.06 0.08 0.1 
六 三 003 上 | et 9 . ts 
0.02- | 改进 MCM 外 7 ” 径 1、3、5 信道 时 域 波形 
0.01 | / Spa a | 
a 3.6 实验 仿真 
9 | 2 六 or3 只 需要 改变 k 的 值 ， 就 可 以 快速 生成 多 个 单 径 自 相关 统计 
ee i 特性 好 ， 且 各 路 径 不 相关 的 瑞 利 衰落 。 采 用 本 文 提出 的 信道 生 
图 5 单 笃信 道 ACF 均 方 误差 a 本 ee 、 
成 方法 , 直接 用 于 MIMO 无 线 信 道 的 仿真 实现 。 仿 真主 要 参数 
3.4 各 k 值 的 MSE 设置 如 表 3 所 示 


mk 


表 3 参数 配置 
收发 天 线 ”调制 方式 ”编码 方式 ”检测 算法 


在 AOA 中 引入 大 值 会 影响 各 径 信道 的 自 相 关 统计 特性 ， 
因此 ， 为 了 更 加 直观 、 全 面 准确 地 了 解 改进 模型 的 统计 特性 ， 
仿真 给 出 了 当 M =5 时 模型 每 一 径 的 自 相关 MSE 的 值 ， 如 医 


所 示 。 BPSK STBC MMSE 
2x2 


CN 


实验 仿真 结果 如 图 8 所 示 。 图 中 给 出 了 改进 MCM 模型 和 

MCM 模型 在 不 同 天 线 数 目下 MIMO 系统 的 误 码 率 。 从 图 中 可 

以 看 出 ， 在 相同 的 信 噪 比 下 ，MIMO 天 线 数目 从 1x1 增 加 到 

2x2 时， 系统 的 误 码 率 依次 下 降 ， 如 实际 相符 ; 在 BER 为 0.1 

s 处 ， 天 线 数目 为 1x1 时 ， 改 进 MCM 模型 相对 于 MCM 模型 有 

2 dB 的 增益 ; 天 线 数目 为 2x2 时 ,改进 MCM 模型 相对 于 MCM 

同 芭 尖 福 这 和 祖坟 守 祷 疾 模型 有 2.1 dB 的 增益 。 在 改进 信道 模型 下 ， 系 统 误 码 率 性 能 得 
到 了 优化 。 图 中 横 坐 标 表示 信 噪 比 ， 纵 坐标 表示 误 码 率 。 
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图 8 各 模型 不 同 天 线 数目 MIMO 系统 误 码 率 
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本 文 提出 了 一 种 基于 SOS 的 多 径 瑞 利 衰落 信道 的 改进 模 
型 ， 在 MCM 的 多 普 勒 频率 中 引入 多 径 因 素 之 后 ， 改 进 了 模型 
的 AOA 取 值 方法 。 在 实现 开销 相同 的 情况 下 , 与 MCM 相 比 ， 
改进 模型 各 径 的 自 相关 函数 的 精度 有 了 极 大 的 提高 ， 且 优 于 其 
他 的 改进 模型 。 改 进 的 模型 可 以 快速 生成 多 个 不 相关 瑞 利 衰落 
言 道 , 可 直接 用 于 MIMO 无 线 信道 和 WSSUS 等 信道 的 仿真 实 
现 ， 拓 展 了 研究 无 线 通信 系统 性 能 时 认识 、 利 用 信 
在 多 普 勒 频率 中 引入 多 径 因 素 提高 各 信道 自 相关 性 能 的 同时 ， 
如 何 进一步 增加 生成 不 相关 波形 的 个 数 ， 并 在 此 信道 下 研究 通 
言 系统 的 性 能 ， 是 下 一 步 研究 的 方向 。 
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